Entanglement Entropy of Coupled Conformal Field Theories and Fermi Liquids 
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In this paper we calculate the entanglement Renyi entropy of two coupled gapless systems in 
general spatial dimension d. The gapless systems can be either conformal field theories (CFT), or 
Fermi liquids. We assume the two systems are coupled uniformly in a ft— dimensional submanifold 
of the space, with < h < d. We will focus on the scaling of the Renyi entropy with the size of 
the system, and its scaling with the inter-system coupling constant g. Three approaches will be 
used for our calculation: (1) exact calculation with ground state wave-functional, (2) perturbative 
calculation with functional path integral, (3) scaling argument. 

PACS numbers: 



I. DEFINE THE PROBLEM 

The entanglement entropy of a system %, refers to the 
entropy of the reduced density matrix of subsystem A, 
after tracing out the rest of the system T-L/A. Usually 
A and H/A are simply separated spatially. If % is a 
conformal field theory (CFT), it is well- understood that 
for dimension d > 1, the leading contribution to the en- 
tanglement entropy Sa is proportional to the size of the 
boundary of subsystem A^; while at d = 1 , Sa is propor- 
tional to ln(Lx), La is the size of Ak~— . If rl is a Fermi 
liquid, then the leading order entanglement entropy is 
Sa ~ L A ~ X \og(LX) with arbitrary dimension g£~— . 

In this work, we study the entanglement entropy be- 
tween two coupled systems A and B in general dimension 
d, namely the entanglement entropy of the reduced den- 
sity matrix pa, after tracing out B. The entire action of 
the system can be schematically written as 

S= I ' drd d xC A + C B + gO A {x)0 B {x)V{x). (1) 



Oa and Ob are normal ordered operators {(O a) — 
(Ob) = 0) of A and B respectively, V(x) is a function of 
space. The entanglement entropy is in general a compli- 
cated function of g (even for most of the exactly solvable 
cases), so we will always assume a weak coupling g, and 
focus on the scaling of the entropy with weak g. Our 
current work will focus on the Renyi entropy, which is 
defined as 



1 - n 



log(tr^) 



(2) 



We take n = 2 as an example hereafter, but the results 
of our paper are insensitive to n. 

Our system is defined on a d— dimensional space with 
coordinates (x\, • • • , Xd), and —L < Xi < L. We assume 
A and B are coupled uniformly in a h— dimensional sub- 
manifold M. of the d dimensional space, i.e. 

V(x!, ■ ■ ■ x h , 0, • • • , 0) = 1, V(x) = otherwise(3) 

When h = d, the coupling is uniform in the entire space, 
while when h = the coupling is only through a point 
contact. 



This paper is outlined as follows: In section II we will 
introduce a general functional path integral and pertur- 
bation theory for the entanglement entropy, and it will 
be applied to the coupled Fermi liquids in section III. 
In section IV, we will study the entanglement entropy 
of coupled CFTs. Besides the perturbation theory, we 
will also introduce a scaling argument to understand the 
qualitative behavior of the entanglement entropy. The 
results with h = d will also be checked with exact calcu- 
lations based on the ground state wave functionals of the 
coupled CFTs. 



II. FUNCTIONAL PATH INTEGRAL AND 
PERTURBATION THEORY 

In most cases, the entanglement entropy between cou- 
pled systems cannot be calculated exactly. In this section 
we will introduce a general perturbation theory based on 
functional path integral formalism for the entanglement 
entropy. 

First of all, the reduced density matrix p A of subsystem 
A at zero temperature is 



hm — 7-57 

/3->oo 7\P) 
AB 



(<Pa\PaWa) = tT B[PAB\ 



x exp 



drd d x C(ipa,<Pb) 



(4) 



Z ao, is the partition function of the entire system: 



Z AB = I D( PaD(pb exp 



drd d x C(ipa,<Pb) ■ (5) 



In the Lagrangian, if ifA and ip& are boson (fermion) 
fields, they are periodic (antiperiodic) in imaginary time 
t € (0,/3). In this section, ifA and tps w iU be taken as 
boson fields for example. 

In order to calculate the Renyi entropy with n = 2, we 
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need to evaluate the following quantity: 
tr[p^] = lim 



AB 



Dtp A Dip B 



f f 2fs 

exp[— / d d x I drC A (if A 



p-e 



20 



dr{C B {^ B )+gO A B V{x))}. (6) 



JO JP+e/ 

In the numerator of this equation, we will keep 
MO) = ¥U(2/3), 



ip B (0) = ip B (/3 - e), ip B (/3 



<PbW)- (7) 



e is an infinitesimal positive number. Field ip B is inte- 
grated out in imaginary time segments r £ (0, j3 — e) and 
t 6 (/3 + e, 2/3) separately. 

The periodicity difference of f A {r) and <^b(t) is the 
key of this calculation. Although we are always consid- 
ering the case with zero temperature, the limit j3 — > +oo 
should be taken after all the calculations with finite (3. 

Now we try to calculate the entanglement entropy with 
perturbation of g. Since we assumed that O a and O b are 
both normal ordered, the first order perturbation of g of 
Eq. [B] vanishes. Expanding Eq. [5] to the second order of 
g, we obtain 

log(tr^)- lim (S 1 -S 2 ), 

P-t+oo 

Si ~ g 2 ^ +J P ^j dT l dT 2 d d x l d d x 2 V{x l )V{x 2 ) 

X (j\j({Tl,Xx, T2,X2)G B g(T 1 ,Xl, T2,X 2 ); 



(8) 



S 2 ~ 2g 2 / dT 1 dT2d d x 1 d d x 2 V(x 1 )V(x 2 ) 
Jo 

X Gaa(ti,Xi, T2,X 2 )G^(ti,Xi, T 2 ,X 2 ). 

G AA an d G BB are correlation functions: 

Gaa(tx,Xi, t 2 ,x 2 ) = (O a (t 1 ,x 1 )O a (t 2) x 2 )), 

G BB (n,xi, t 2 ,x 2 ) = (O b (t 1 ,xi)O b (t2,x 2 )), 

,Xl, T 2 ,X2 

,x%, t 2 ,x 2 ). (9) 

Notice that there are two different periodicities in these 
correlation functions. If A and B are both bosonic op- 
erators, then in the frequency space, G^> has Matsubara 
frequency 2irm/j3, while G^ 2 ^ has frequency 2-7rm/(2/3). 
For example, if A and B are CFTs with z = 1, and we as- 
sume operator A ( B ) has scaling dimension A_4(g), then 
G^ reads 



G 



iP) 

AA 



(&B)(0>0> T i x ) 



E 



PL d ^ \oj 2 + k 2 

uj.k 



i(d+l-2A A(B) ) 



iu)T-\-ik-x 



(10) 



Eq. |B] and Eq. [5] were formulated for Renyi entropy 
with n — 2 only, but their generalization to arbitrary n 
is straightforward. 



III. ENTANGLEMENT ENTROPY OF 
COUPLED FERMI LIQUIDS 

A. Uniform tunnelling 

In this section we will consider the entanglement en- 
tropy of coupled Fermi liquids. The simplest situation 
that we can start with, is that A and B are free Fermi 
gases with S z = ±1/2 respectively, and they are cou- 
pled together through a uniform transverse magnetic field 
HS X . Our goal is calculate the entanglement entropy of 
S z — 1/2 fermions, after tracing out the S z = —1/2 
fcrmions. With uniform magnetic field, this system can 
be trivially solved, and the reduced density matrix p A 
is a simple direct product of the density matrix at each 
momentum k: 



PA 



k 



(11) 



If both S x = ±1/2 spin states are occupied or unoccu- 
pied, PA,k is a pure state density matrix. While if only 
one of the spin states is occupied, p A _^ is maximally en- 
tangled: 



1 f 

pam = -4 >t 



Ck,, + \ 



(12) 



Therefore only the states with energy £/ — H/2 < e < 
Ef + H/2 contribute to the entanglement entropy. Hence 
the entanglement entropy should scale as 

S A ~Af(e f )\H\L d . (13) 

N(sf) is the density of states at the Fermi surface. 

B. Point contact tunnelling 

Now suppose spin up and down fermions are coupled 
through a static polarized magnetic impurity at r = 0: 
HS X (Q), this impurity tunnels spin up and down fermions 
through the point contact at r = 0. The perturbation 
formalism developed in the previous section is applicable 
here, as long as in the calculation we keep the Matsubara 
frequency for spin up (A) and down (B) fermions as 



U) A 



7r(2m + l) 



L) B 



?r(2n+ 1) 



(14) 
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Notice that the difference between Matsubara frequencies 
uj a and uj b leads to 



(-1)' 



dr exp (i(uj A — wb)t) = 
u u A - U B 



(15) 



which contrasts the delta function in the usual case. 

The leading order contribution to the entanglement en- 
tropy is a straightforward application of Eq. [8j and it 
leads to the following results: 

log (trp A ) = ptofJS! - S' 2 ), 



E E 



2 L- 2d 



1 



^(UU-WB) 2 
1 



i^A - £k,A + £f i^B - Zk.B + £/ 



2^ 2 E E ^ 

1 



2 c? 



1 



£fe,yt +£/ KjJ — £fc,B + £/ 



(16) 



Frequency u takes the usual values 7t(2to + l)//3. 

Correct evaluation of the frequency and momentum 
summation in Eq. 1161 leads to the following result: 



f + OO r-Zf r-Zf 

S A ~ H 2 ( / de A / de B + / de A 
Jo Jo 

jN(e A )N(e B ). 



dsr 



{£A - SB, 



(17) 



Since the density of states is a constant close to the Fermi 
surface, this integral is logarithmically divergent when 
e A and eq are close to Fermi energy £/. This logarithmic 
divergence will be cut-off by 1/L, thus the final result of 
the entanglement entropy is 



S A ~H 2 (A/-( £/ )) 2 log(L). 



(18) 



In one dimension, the Fermi liquid becomes Luttinger 
liquid, which is a CFT. Uniform magnetic field and point 
contact single fermion tunnelling have scaling dimensions 
A = 1 and respectively on a free fermion Luttinger liq- 
uid CFT, i.e. the point contact single fermion tunnelling 
is a marginal perturbation on the free fermion Luttinger 
liquid. Later we will see that the results in Eq. Q2] and 
Eq. [TS] are consistent with our general results about CFT 
with d = 1. 

Many aspects of the Fermi liquid theory can be viewed 
as infinite number of one dimensional fermions moving 
along the radial direction, thus it is not surprising that 
the entanglement entropy of Fermi liquid at higher di- 
mension is qualitatively equivalent to one dimensional 
free fermions. The connection between the Fermi liquid 
and one dimensional CFT was also used to understand 
the ordinary entanglement entropy of Fermi liquid 7 . 



C. h dimensional tunnelling 

Now Let us assume the S z = ±1/2 fermions 
are coupled through a transverse magnetic field on a 
h— dimensional submanifold Ai of the space (Eq. [3]). If 
we take the simplest quadratic fermion dispersion, the 
second order perturbation in Eq. [S] gives the following 
result: 

S A = Si — S'2, 



Si - H 2 £ £ 

Wi.WB ki,kj^j ,ks,j 
1 



2 L 



2h-2d 



p 2 {w A -u: B y 



1 



s> ~h 2 y, E L 

u k i: k A j ,k Bd 

1 



2h-2d 



- Eti k i - Y? =h+l k A,j + £ f 

1 



IU1 - 



Ei=i k f - J2j=h+i k B,j + £ f 



(19) 



When h < d, this integral is always logarithmically 
divergent, thus the logarithmic contribution persists (at 
least to the second order perturbation) as long as h < d: 



S A ~ H 2 L h log(L) 



(20) 



IV. ENTANGLEMENT ENTROPY OF 
COUPLED CONFORMAL FIELD THEORIES 

If the two coupled systems A and B are both CFTs, 
the entanglement entropy due to coupling g will obvi- 
ously depend on the scaling dimension A of the cou- 
pling constant g. If the scaling dimensions of A and 
O b are A_4 and Ag respectively, then the dimension of g 
is A = h + z — A A — A/3, z is the dynamical exponent. 



A. Exact calculation with ground state 
wave-functionals 

We will first consider the following theory 



C = 



\d T <P A t\ 2 + \d T V B i\ 2 + \k z <P A t\ 2 + \k z V B>% \ 



+ ^n^a + ^gl 2 +9bh™\v A$ -V BS \ 2 . (21) 

Both A and B are free boson theories, and the coupling 
between them is uniform in the d— dimensional space. 
Here a and b are both dimensionless constants, and q is a 
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small coupling constant. By adjusting number to, we can 
tune the scaling dimension of g: A = 2z — to. Since the 
entire Lagrangian Eq. [21] is quadratic, the entanglement 
entropy Sa can be calculated exactly. 

The entanglement entropy Sa can be calculated in the 
same formalism as Rcf. 8 , where a marginal coupling be- 
tween two Luttinger liquids was considered (z = 1, to = 
2). Since the exact ground state wavefunctional of field 
ifA and ips can be written down exactly, one can directly 
calculate the entropy with the exact reduced density ma- 
trix. The coupling in Eq. is uniform in space, it only 
couples tpA and Lps modes with the same momentum, 
hence the entropy is a simple sum of the entropy of cou- 
pled harmonic oscillators at each momentum k. Since 
the exact result is in general a complicated function of 
g, we will focus on the leading term after expanding the 
exact result with small g. 

The leading order results with different choices of A 
are summarized as follows: 



(1), A = d/2, 



1 



S A ~g 2 \og( — )L d ; 

\9\ 

(2) , A > d/2, 

S A ~9 d/A L d , 

(3) , A < d/2, 

S A -9 2 L d . 



(22) 



Notice that d/2 is the critical value of A, when A < d/2 
the leading term of the Renyi entropy is always g 2 L d , 
whether g is relevant or not. When A = d/2, the lead- 
ing order entropy acquires a logarithmic correction. In 
Eq. [22j we have assumed that the infrared cut-off of the 
system is g instead of 1/L, i.e. L > y _1 / A . If L < 3 _1 ^ A , 
the argument of the logarithmic function in Eq. [2D is re- 
placed by L. 

It is known that the subleading correction to the 
boundary law of the ordinary entanglement entropy con- 
tains important information about the CFT— In fact, 
in Eq. [22j in addition to the leading term proportional 
to the volume of the system, there are also subleading 
terms. The subleading terms can be calculated conve- 
niently for free boson theory Eq. [5T] at d = 1, and we 
summarize our results here: 



subleading 



log 



'(VE + Vb) 4 

lQab 



, A>0, 



blcading 1 



blcading ' 



Alog(L), 6 = 0, A>0, 
g 2 , A = 0. 



(23) 



When A < 0, there is no subleading term at order 0(L°). 

Interestingly, when the coupling is relevant i.e. A > 0, 
the most generic subleading contribution to the entropy 



approaches a constant when g — > 0, as long as we take the 
limit L — > oo first. If b — (i.e. the coupling only affects 
mode ip A + ¥B, while ip A - <p B is still gapless), 5 sub icadin g 
is logarithmic of the system size, and its coefficient is a 
universal constant proportional to the scaling dimension 
A, but it is independent of the magnitude of g. 

The logarithmic subleading term in Eq. [231 may have 
generalizations to other CFTs in one dimension. The 
universal coefficient of the logarithmic term might be re- 
lated to the central charge of the CFT. We will study 
this general theory in future. 



B. Perturbative calculation 

We can also apply our perturbative formalism Eq. [5] 
to the coupled CFTs. In the momentum and frequency 
space, the entropy is evaluated as 



log(tr^) = lim {S' x -S' 2 ), 



si 



' E E 



4 L 



2h-2d 



1 



1 

l w B + Z^=l fe i + 2^j=h+l K B,j> 2 



- .<r'E E 1 



j=h+l -J2,jj 
2h-2d 

1 



i 



(- 2 + Eli kf + E?=m-i fci ;i )^ +1 - 2Aa > 



(24) 



We have taken z = 1 as example. A-a(b) is the scaling 
dimension of Oa(b)- I n this equation, 

2tt(to+1/2) 2im , , 

f>A = V p ' - , w, u B - — . (25) 

Correct evaluation of the summation in Eq. [24] will lead 
to the results consistent with the exact results Eq. [521 In 
this calculation, one should always take the limit f3 — > 
oo before the limit L — > oo. When h = d (two CFTs 
are coupled uniformly in the entire space), the leading 
contribution to Eq. [M] is 



s A ~g 2 L d 



l 



^,2rf+2-2A^-2Ai 



-d d k. 



(26) 



The scaling dimension of g is A = d + 1 — A a — Ag. 
If we take A = d/2, this integral gains a logarithmic 
contribution. Since A > 0, the higher order perturba- 
tion will acquire stronger and stronger infrared diver- 
gence. If L < g~ x ^ A , this logarithmic divergence is 
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cut-off by l/L; if L > .g _1 / A , we expect the summa- 
tion of the perturbation series will eventually be cut- 
off by length scale <? _1 / A , so the final answer should be 
Sa ~ 9 2 l°g(l/l5l)^ d ! which is consistent with Eg. |2"!?1 

Similarly, if we take < h < d, then when A = h/2 
the perturbation theory gives the logarithmic term 

^~5 2 log(^|)^, (when A = h/2). (27) 

Now suppose we take h — (point contact), the scaling 
dimension of g is A = 1 — A A — Ag. For simplicity, we 
assume that A A = A B = |(1 - A). Then Eq. [2j] is 
evaluated as 

- 9 2 f j^dk. (28) 

When A < 0, the point contact is irrelevant, and this 
integral gives a constant result S A ~ g 2 , with a nonuni- 
versal coefficient. When A = (marginal point contact), 
we obtain a logarithmic contribution, regardless of the 
total dimension d: 

S A ~g 2 log(L). (29) 

In this case, since the point contact is marginal, the 
higher order perturbation has the same logarithmic di- 
vergence as the leading order, hence we expect this loga- 
rithmic divergence persists even after we sum the entire 
series, and it can only be cut-off by L. 



C. Scaling Argument 

In this section we will try to understand the results 
obtained in the previous two subsections using a sim- 
ple scaling argument. The argument in this section is a 
generalization of RefJ^, where a scaling argument was 
introduced to explain the logarithmic contribution of the 
ordinary entanglement entropy of Id CFT. Scaling ar- 
gument was also introduced to understand the entangle- 
ment entropy close to finite temperature critical points^. 

We first note that, if the two coupled systems are 
gapped, the Renyi entropy is obviously proportional to 
the volume of the coupled submanifold, and it should 
scale as g 2 for weak coupling: S A ~ g 2 L h . This is be- 
cause S A is positive definite, hence the leading contribu- 
tion should be an even function of g. Also, the entan- 
glement entropy vanishes in the decoupling limit g — > 0. 
Thus the leading analytic contribution from the coupling 
should be g 2 . 

Now we consider two CFTs coupled on a 
h— dimensional submanifold A4 of the space, and 
the entanglement entropy is collected while coarse- 
graining the system. At length scale I, the size of the 



coupled subsystem is effectively L h /l h . Within length 
scale interval cZ(logZ), the entanglement entropy is 

d(S A )~^g 2 d(logZ). (30) 

gi is the effective coupling constant at length scale I, and 
as long as gi is weak, gi ~ gl A . Thus if A > 0, the total 
entropy is 

r l=Min[\g\- 1/& , L] jh 

Sa - / d (log I)— g 2 l 2A . (31) 

Jl=a 1 

This integral gives us the following results for general 
d> h> 0, with different choices of A: 

(1) , A = h/2, 

S A ~g 2 log(±-)L h , Lyyg- 1 ^; 

m 

S A ~g 2 L h log(L), L«g- 1/A ; 

(2) , A > h/2, 

S A ~g h ' A L h , L»g-^ A - 
S A ~g 2 L 2A , L« 5 - 1 /A. 

(3) , A < h/2, 

S A - g 2 L h . (32) 

Now the critical value of A becomes h/2. When h = d, 
and L > g~ 1/A , Eq. [521 reduces to the results Eq. |2"2I 
obtained from exact calculations. 

In the special case with h — 0, i.e. the CFTs are cou- 
pled through a point contact, the integral in Eq . [3T1 gives 
us the following results: 

A = 0, S4~.g 2 log(L), 

A>0, S A ~C, L».g- 1 / A , 

S A ~g 2 L 2A , L«g-^ A , 

A<0, S A ~g 2 . (33) 

C is a constant which does not scale with g. When A = 0, 
namely the case with a marginal point contact, the lead- 
ing contribution to the entanglement entropy is a loga- 
rithmic term, which is independent of the spatial dimen- 
sion d. This conclusion confirms our calculation in the 
previous subsection (Eq. [551 Eq. |2"9"|) . and confirms the 
calculation for Fermi liquid with a point contact (Eq. [18]) . 

This simple scaling argument should be precise if g is 
irrelevant or marginal, for arbitrary d and h. If g is rel- 
evant, the integral of length scale in Eq. [31] was taken 
only from a (lattice constant) to the scale where g be- 
comes nonperturbative. It seems like we have ignored 
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the entropy contribution after g becomes nonperturba- 
tive. To understand this problem, we need to know the 
long wavelength properties of the system when g is rele- 
vant, and there are two possibilities: 

1. In most cases, a relevant coupling g opens up a 
gap (or local gap) for CFTs A and B, namely all the 
correlation functions G(t± — T2, x\ — x%) with X\,x% G M. 
decays exponentially when |n — rg| — > 00. In this case, 
the system can be driven into either a direct product state 
between A and 6, or a maximally or partially entangled 
state between A and B with a saturated entanglement. 
For instance, if A and B are free boson fields, and in the 
submanifold M. they are coupled as —ag(ipA + Ve) 2 — 
bg{ipA ~ <Pb) 2 , then as long as ag > and bg > 0, the 
relevant coupling g will keep tpA = ipjs = in M, hence 
the system becomes a direct product state between A 
and B in the long wavelength limit. Then the entropy 
with length scale I > g~ x ^ A is ignorable, and there is no 
correction to Eq. [3T] 

When g drives the system into a maximally or partially 
entangled state, then we need to add another contribu- 
tion to the entropy, which is 



5" 



A 



{ g-l/A )h 



= g h / A L>> 



(34) 



I and 



With this extra contribution, our results in Eq. 
Eq. [3j still hold. 

2. If some of the correlation functions in the coupled 
submanifold M remain power-law even with relevant g, 
then there is a residual scaling invariance after g becomes 
nonperturbative. In this case, we need to include the 
following extra contribution to Eq. [3T] 



S' 



A 



l=L 



-l/A 



dOogO-pr 



(35) 



This integral does not modify any of the leading order 
terms in Eq. [32] with h > 0, but it leads to a logarithmic 
contribution to Eq. [33] with A > and h = 0, i.e. it only 
affects the case with a relevant point contact coupling 



between A and B. For instance, in Eq. [21] although the 
mode ip ^ — tpB remains gapless when 6 = 0, the exact 
results Eq. [22] always agree with Eq. [32] obtained from 
scaling integral Eq. [3T]for h = d, no matter b = or not. 



V. SUMMARIES AND EXTENSIONS 

In this work, we studied the entanglement entropy of 
coupled Fermi liquids and CFTs. Three different meth- 
ods were used for the calculation: perturbation the- 
ory, scaling argument, and exact ground state wave- 
functional. These three approaches are consistent with 
each other for all the cases that we can check. 

It has been demonstrated that the holographic method 
is a very powerful way of calculating the entanglement 
cntropy2ii£ of CFT, assuming there is a bulk AdS space 
duality of the boundary CFT. The ordinary entangle- 
ment entropy is related to the area of the minimal surface 
of the bulk AdS space. In future, we will try to develop 
a holographic formalism to produce the results in the 
current paper. Since AdS/CFT duality effectively "ge- 
ometrizes" the RG flow at the boundary CFT theory, we 
expect the holographic calculation of the entanglement 
entropy to be qualitatively equivalent to the scaling ar- 
gument discussed in this paper. 

Significant progresses have been made in numeri- 
cal simulation of quantum many-body states. For in- 
stance, the multi-scale entanglement renormalization 
ansatz (MERA) is especially powerful in simulating one 
dimensional CF T 14 ' 15 . In future, it will also be interest- 
ing to verify the conclusions in our paper numerically. 
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